Time scale competition leading to fragmentation and recombination transitions 

in the co-evolution of network and states. 
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We study the co-evolution of network structure and node states in a model of multiple state 
interacting agents. The system displays two transitions, network recombination and fragmentation, 
governed by time scales that emerge from the dynamics. The recombination transition separates a 
frozen configuration, composed by disconnected network components whose agents share the same 
state, from an active configuration, with a fraction of links that are continuously being rewired. 
The nature of this transition is explained analytically as the maximum of a characteristic time. 
The fragmentation transition, that appears between two absorbing frozen phases, is an anomalous 
order-disorder transition, governed by a crossover between the time scales that control the structure 
and state dynamics. 
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I. INTRODUCTION 

Recent findings in the topological characterization of 
real complex networks have triggered a theoretical un- 
derstanding of diverse complex systems. A great deal 
of effort has been devoted to the modelling of complex 
networks from a topological point of view, and to the dy- 
namics on different classes of networks In the latter, 
the evolution of the states of the nodes is usually assumed 
much faster that the time characterizing the network dy- 
namics. Less effort has been devoted to the understand- 
ing of the entangled co-evolution of network structure 
and state dynamics, i.e., how the structure of a network 
affects the dynamics on it and vice versa. For instance 
in a social system, individuals shape their opinions de- 
pending on their neighbors' opinions, but simultaneously, 
the individuals' opinions affect with whom they interact 
@> B BL IS B 0] ■ The network of interactions co-evolve 
with the dynamics of opinions. 

The influence of network topology has been analyzed 
extensively for several models of consensus formation 
[1)3- A main question is related to the mechanisms 
and network topologies that lead to consensus, i.e., all 
agents having the same state. Based on social pressure, 
how state changes, and homophily, the tendency of indi- 
viduals to interact with similar others, the Axelrod model 
[Io| represents a paradigmatic model displaying an order- 
disorder noncquilibrium transition [111 ]. When the net- 
work of interaction is regular, random or small-world, the 
system orders if the degree of initial disorder is below a 
critical value [ll], [IH, 0, [Hj]. In the ordered phase, a 
domain (set of connected agents with the same state) of 
the order of the system size spans the network, while in 
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the disordered phase, many small domains are formed. 
Motivated by the mechanisms of social pressure and ho- 
mophily, in this paper we analyze the co-evolution in the 
Axelrod model between the network of interactions and 
the state dynamics of the nodes, by allowing interaction 
links to be rewired depending on the state of the agents 
at their ends. 

II. THE MODEL 

A population of TV agents are located at the nodes 
of a network. The state of an agent i is represented 
by an F-component vector c^/, / = 1,2, . . . , F, and 
i = l,2,...,N, where each component represents an 
agent's attribute. There are q different choices or traits 
per feature, labelled with an integer Uif £ {0, . . . , q — 1}, 
giving raise to q F possible different states. 

Initially agents take one of the q F states at random. 
In a time step an agent i and one of its neighbors j are 
randomly chosen: 

1. If the agents share m > features, they interact 
with probability equal to the fraction of shared fea- 
tures, i.e., the overlap (m/F). In case of interac- 
tion, an unshared feature is selected at random and 
i copies j's value for this feature. 

2. If the agents do not share any feature, then i dis- 
connects its link to j and connects it to a randomly 
chosen agent that i is not already connected to. 

Step 1 describes the original Axelrod dynamics: Alike 
agents become even more similar as they interact, in- 
creasing the probability of future interaction. Step 2 im- 
plements the network co-evolution: incompatible agents, 
i.e., agents with no features in common, tend to get dis- 
connected. We have performed extensive numerical sim- 
ulations to study the behavior of the system as the con- 
trol parameter q is varied, for different population sizes 
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FIG. 1: Average relative size of the largest network compo- 
nent (circles) and largest domain (solid line) in the stationary 
configuration vs q, for N = 2500, averaged over 400 realiza- 
tions. The vertical lines at q c = 85 and q* = 1875 indicate the 
transition points between the different phases. Inset: fluctu- 
ations are maximum at the critical point q c . 



N, number of features F = 3, and starting from a ran- 
dom network with average degree (k) = 4. The results do 
not depend on the initial network topology, because the 
repeated rewiring dynamics leads to a random network 
with a Poisson degree distribution. 

The model displays two transitions, both very different 
in nature. The first one is an order-disorder transition 
at q = q c between two frozen phases, associated with the 
fragmentation of the network. The dynamics leads to the 
formation of network components, where a component is 
a set of connected nodes. In a frozen configuration agents 
that belong to the same component have the same state. 
For q < q c (ordered phase I), the system reaches a con- 
figuration composed by a giant component of the order 
of the system size, and a set of small components; for 
q > q c (disordered phase II) the large component disin- 
tegrates in many small disconnected components. The 
second transition, related with network recombination, 
occurs at q = q* between the phase II and an active 
phase III, where the system reaches a dynamic configu- 
ration with links that are permanently rewired. 



III. FROZEN PHASES (I AND II) 

For small values of q in phase I, the average size of 
the largest network component S in the final configura- 
tion is of the order of the system size N (Fig. |T|) , due to 
the high initial overlap between the states of neighboring 
agents. As q increases inside phase I, the initial overlap 
decreases and S also slowly decreases. For larger values of 
q (phase II), many distinct domains are formed initially 
inside the components which break into many small dis- 
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FIG. 2: Size distribution of network components for TV = 
2500 and values of q (a,b) below, (c) at, and (d) above the 
transition point q c ~ 85. The dashed line represents a power 
law with exponent —1.3 ± 0.02. 



connected components and, as a result, S reaches a value 
much smaller than N (network fragmentation). During 
the evolution of the system, a network component can 
have more than one domain. However, in the final con- 
figuration of phases I and II, one component corresponds 
to one domain. 

The transition point from phase I to phase II is de- 
fined by the value q = q c for which the fluctuations in 
S reach a maximum value. This value corresponds to 
the point where the order parameter S suffers a sudden 
drop. For N = 2500 we find q c = 85 ± 2 (see Fig. [TJ. 
To further investigate this transition point we calculated 
the size distribution of network components P(s) (see 
Fig. (2). For small q, P(s) shows a peak that corresponds 
to the average size of the largest component S, and has 
an exponential decay corresponding to the distribution of 
small disconnected components. The peak at S decreases 
as q increases, giving raise to a power law decay of P(s) 
at q c ~ 85, a signature of a transition point [l3l [l5j. 

The behavior of the order parameter S/N depends on 
system size (Fig. [3]). When both axis are rescaled by 
N~ a with a = 0.82 ±0.01 the data collapses for q smaller 
than q c . This implies that q c increases with the system 
size N as q c ~ A^", and suggests a scaling relation for 
S in the ordered phase S = N a f(N~ a q), where /(■) 
is a scaling function. The scaling relation implies that 
the discontinuity disappears in the large N limit. Thus, 
not only the transition point diverges as q c ~ N 82 but 
also the amplitude of the order parameter S/N ~ N~°- ls 
vanishes as N goes to infinity. 

These results, maximum of the fluctuations, data col- 
lapse and distribution of component sizes, identify q c as 
the critical point of the transition. The fact that the ex- 
ponent of the size distribution is smaller than 2 and that 
the discontinuity of the order parameter tends to as 
system size increases, suggest that in the large A^ limit 
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FIG. 3: Average relative size of the largest network compo- 
nent S/N vs q for system sizes N = 2500, 6400, 10000 and 
14400 (left to right), averaged over 400 realizations. Inset: 
Finite size scaling. The data collapses below the scaled tran- 
sition point Q c = q c N~ a , with a — 0.82. 



FIG. 4: Time evolution of the density of domains rid (solid 
line) and network components n c (dashed line), for N = 2500 
and values of q (a) below, (b) at, and (c) above the transition 
point q T — 194 ± 10. Each inset is a zoom of the region that 
shows the approach to the final configuration. 



the transition becomes continuous with q c — > oc. 



IV. ACTIVE PHASE (III) 

We analyze this phase by looking at the rewiring dy- 
namics. A link that connects a pair of incompatible 
agents is randomly rewired until it connects two com- 
patible agents, i.e., agents with at least one feature in 
common. If the number of pairs of compatible agents 
L c is larger than the total number of links in the system 
(k)N/2, this rewiring process continues until all links con- 
nect compatible agents. Later, the system evolves until 
each component constitutes a single domain, the frozen 
configurations reached in phases I and II. If, on the con- 
trary, L c gets smaller than (k)N/2, the system evolves 
connecting first all L c pairs by links. The state dynam- 
ics stops when no further change of state is possible, but 
the rewiring dynamics continues for ever with approxi- 
mately (k)N/2 — L c links that repeatedly fail to attach 
compatible agents. This is the active configuration ob- 
served in phase III. Thus, in contrast to phases I and 
II, in the stationary configurations of phase III there are 
typically more than one domain per component. The size 
of the largest component abruptly increases at q* indicat- 
ing that a giant component reappears (network recombi- 
nation), while the size of the largest domain continues 
decreasing (see Fig. [TJ. 

To estimate q* , we will assume L c constant during the 
evolution as for q large the state of the agents does not 
evolve much. Thus, with the ansatz L c ~ L c {t — 0) ~ 
N{N - 1) [1 - (1 - l/q) F ] ~ N 2 F/2q, for q > F and 
N ^> 1, the condition at the transition point is 



N 2 F/2q* ~ (k)N/2, so that 



is the value of the recombination transition point. For 
the values considered here we obtain q* — 1875 in good 
agreement with the numerical results (Fig. [1]). 

V. DYNAMIC TIME SCALES 

In order to understand the final structure of the net- 
work in phases I and II we analyzed the time evolution 
of the nodes' states and interaction links. In Fig. 0] we 
plot the time evolution of the density of network com- 
ponents, n c — number of components/A^, and domains, 
rid =number of domains/A^, averaged over 1000 realiza- 
tions, and for three values of q. An interesting quantity 
is the average time to reach the final frozen configuration 
t. If Td [t c ] is the average time at which rid [n c ] reaches 
its stationary value, then r is largest between Td and t c 
(see Figs. 0] and [5]). The curves for Td and t c as a function 
of q cross at a value q T (Fig. [5]) . As we shall see, q T iden- 
tifies a transition between two different dynamic regimes 
for the formation of domains and network components, 
that lead to the frozen configurations observed in phases 
I and II (Figs.i] and EJ. 

For q < q T , the dynamics causes the network to break 
into a giant component and small components. Due to 
the initial overlap between the states of the agents in- 
side each component, the network stops evolving at a 
time t c where n c reaches its stationary maximum value 
(see Fig. 2^,). After this stage, domains compete inside 
each component, until only one domain occupies each 
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FIG. 5: Convergence times r c (circles), (squares) and r 
(solid line) vs q for N = 2500 and averaged over 500 config- 
urations. The result from Eq. ((5} (dashed line) is compared 
with t and r c for q > q T = 194 ± 10. Inset: relative differ- 
ence between the transition points q T and q c vs 1 /N in log-log 
scale. The dashed line has slope 0.30 ± 0.01. 



component. The approach to the frozen configuration 
is controlled by the coarsening process inside the largest 
component, whose structure is similar to a random net- 
work due to the random rewiring dynamics. Then, given 
that the dynamics of the last and longest stage before 
reaching consensus inside this component is governed by 
interfacial noise as in the voter model [l(|, r is expected 
to scale as the size of the largest component r ~ S [171 ]. 

For q > q T , there is a transient during which rid de- 
creases, indicating that, in average, domains grow in size 
(see Fig. Hfc). At time Td, rid reaches a stationary value 
when the overlap between distinct domains is zero. At 
this stage domains are still interconnected by links be- 
tween incompatible agents. As domains progressively 
disconnect from each other n c increases. When finally 
the links connecting incompatible agents disappear all 
domains get fragmented, n c equals rid and the system 
reaches its final configuration in a time r = t c . 

At q T , the time scale governing the state dynamics is 
the same as the time scale governing the network dy- 
namics Td — t c (Figs. |4)d[5]). Thus, it indicates that the 
ordered phase I is dominated by a slower state dynamics 
on a network that freezes on a fast time scale, while in 
the disordered phase II the state dynamics freezes before 
the network reaches a final frozen configuration. Even 
though q T is found to be larger than the critical point 
q c , the relative difference between q T and q c decreases 
with N (inset of Fig. [5]), suggesting that both transition 
points become equivalent in the large N limit. Thus, the 
competition between the time scales Td and r c governs 
the fragmentation transition at q c . 

In the remainder, we derive an approximate expression 
for r in phase II by studying the decay of the number of 
links between incompatible agents Nq. We shall see that 



this approach unveils the transition from the frozen to 
the active phase, leading to the transition point q*. 

In the continuum time limit, and neglecting the cre- 
ation of incompatible links, No decays according to the 
equation: 

dN 1 2N N c 2N N c 



dt 



-L (k)N N 



(k)N 



(2) 



In a time step At = i, an incompatible link is 
chosen with probability . One of its ends j is moved 
to a random node k. The probability that k is compatible 
to i is j£ , where N c is the number of compatible agents 
to i but still not connected to i. In a mean-field spirit, 
every node has (k) edges that need to be connected to (k) 
different compatible agents. We approximate N c as the 
average number of compatible agents per agent. When 
i attaches an edge to a compatible agent, No is reduced 
by one while N c is reduced by 2/N given that both i 
and k loose a compatible partner. Assuming that the 
set of compatible agents to i remains the same we write 
N c ~ 2^0. _|_ A, where A = N c (0) ~ H^°2 ~ ^£ - (k) for 
N ^> 1 and q F. Substituting this last expression for 
N c into Eq. ([2|) and rewriting it in terms of the transition 
point q* — NF/(k) we obtain 

dN Q _2No_ /2Ao (k) 
dt ~ (k)N \ N q 

Equation ^ has two stationary solutions. For q < q* , 
the steady configuration is Nq — 0, corresponding to the 
frozen phases I and II, while for q > q* the stationary 
solution Nq — 2^-{g. — q*) corresponds to the active 
phase III. We recover our previous result that for q > q* 
the system reaches a stationary configuration with a con- 
stant fraction of links Nq larger than zero, that are per- 
manently rewired. Therefore, the system never freezes. 
Note that in the limit of very large q, all agents are ini- 
tially incompatible, consequently Nq approaches to the 
total number of links ^(k)N. 

Integrating Eq. ([3]) by a partial fraction expansion 
gives 



(3) 



qN 



2(q* - q) 



In 



(4) 



For q < q* , the system freezes at a time r at which 
N ~ 1, thus 



qN 



2(q* - q) 



In 



q + Uk)N(q*-q) 



1 



for q T < q < q* 



(5) 

This result is in agreement with the numerical solution 

(Fig. ED . 

For q > q* , the system reaches a stationary configu- 
ration. Thus we define r as the time at which A^o — 
^(<Z -<?*) + !, then 



2q 



qN 



2(q* - q) 



In 



2q* 



q*({k)N(q 



2q) 



(0) 
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for q* < q < . r decreases with q and it vanishes for 
q > N 2 F/2 where initially all pairs of nodes are incom- 
patible, thus the system starts from a stationary config- 
uration, giving r = for q > ^j^. 

From Eqs. (O and ([6]) we obtain that r reaches a 
maximum value equal to j(k)N 2 at q = q*\ an indication 
of the transition. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the Axelrod model with co-evolution 
of the interaction network and state dynamics of the 
agents. The interplay between structure and dynamics 
gives rise to two different transitions. First, a recombi- 
nation transition between a frozen and an active phase. 
The characteristic time to reach a stationary configura- 
tion shows a maximum at the transition point between 
these two phases. Second, an order-disorder transition 
associated with network fragmentation that appears at 
a critical value q c where the component size distribution 
follows a power-law. Finite size scaling analysis suggests 
that in the large N limit this transition becomes con- 
tinuous with q c going to infinity. The fragmentation is 
shown to be a consequence of the competition between 
two coupled mechanisms, network formation and state 
formation. These mechanisms are governed by two inter- 



nal time scales, r c and Td respectively, which are not con- 
trolled by external parameters, but they emerge from the 
dynamics. For q < q c the network components, that are 
formed first, control the formation of states. For q > q c 
the fast formation of domains shape the final structure 
of the network. 

An important aspect of the co-evolution is that the 
network evolution is coupled to the state of the agents, 
in contrast to other models where the links are severed, 
reconnected and/or appear as a random process inde- 
pendent of the state of the nodes. The robustness of 
the fragmentation and recombination in the presence of 
noise should also be considered [l8j]. Our results provide 
a simple mechanism, co-evolution, that could explain the 
community structure found in general in the analysis of 
complex networks, and in particular of social systems 
[l9| where this mechanism can be understood in terms of 
network homophily [l8| . 
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